Some properties of membranes in nematic solvents 
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The fluctuation spectrum of membranes in nematic solvents is altered by the boundary condition 
imposed on the bulk nematic director by the curved membrane. We discuss some properties of 
single and multi-membrane systems in nematic solvents, primarily based on the Berreman-de Gennes 
model. We show that: membranes in nematic solvents are more rigid and less rough than in their 
isotropic counterparts; have a different Helfrich steric stabilization energy, proportional to d"^, and 
hence a different compression modulus in the lamellar state; and can exhibit phase separation via 
unbinding during a quench into the nematic state. We also discuss the preparation and possible 
experimental effects of nematic-mediated surfactant membrane system. 
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There are many instances in physics where the configuration of a field on a boundary is infiuenced by the fluctuations 
of a conjugate field in the adjacent volume. Examples of this are the Casimir effect, where the confinement of the 
electromagnetic field between conducting plates results in a net attraction jl| : related effects in soft condensed matter 
include a fiuctuation-enhanced interaction between inclusions on a membrane Q and the the interaction between 
surfaces dipped in a structured fiuid due to the change in fluctuation spectrum of the fluid 1^,1). While these are 
entropic in origin, there is another class of fluctuation enhancements due to the energy of deforming the bulk field 
coupled to the boundary field by some anchoring condition. Examples here include: interaction between membrane 
inclusions due to the strain induced in the membrane |^]; the non-analytic contribution to the wetting contact line 
elasticity due to deformations of the adjacent fluid-air interface and the example which we explore here, a non- 
analytic contribution to the free energy of a surface in contact with a liquid crystalline solvent ||^,^ . 

In this article we explore some of the properties of fiuid membranes in contact with nematic liquid crystalline 
solvents, discussing both entropic and energetic effects. By 'membrane' we envision surfactant bilayers arranged in 
the archetypical structures found in surfactant systems: in this work we focus on lamellar phases. It is well known 
that surfactants induce varying degrees and strengths of boundary conditions on the nematic director . Here we 
consider the simple natural case when the mesogenic molecules are strongly anchored by the hydrophobic tails of 
surfactant in the direction along the membrane normal. 

Our starting point is the well-known Berreman-de Gennes model 0,^, which was introduced to describe the 
anchoring energy of a liquid crystal on grooved substrates. Here we consider the 'grooved substrate' to be thermal 
undulations of a bilayer surface, and hence consider some of the consequences of the deformation of the nematic 
director field, induced by these undulations, Fig.l. The principle is this: the equilibrium thermodynamics of a 
membrane-liquid-crystalline-solvent system includes as fluctuating variables both the solvent and the membrane. If 
we are interested in properties of the membrane we can 'integrate out' the solvent degrees of freedom to find an 
effective theory for the membrane thermodynamics. To perform such an integration and speak of a renormalized 
theory of membranes requires a separation of timescales. That is, if we are interested in dynamical properties of the 
membrane, we must ensure that the solvent fluctuations (in this case, the establishment of a deformed director field 
in response to a surface undulation) are much faster than the characteristic membrane decay time. However, if we 
are interested in equilibrium effects, such as fluctuation spectra as would be measured in experiments lasting 'long' 
times, then this procedure is valid. It is with these kinds of experiments in mind that we proceed. Let us first recall 
the fundamental ideas of the Berreman-de Gennes model. 

Nematic Energy — Consider a surface with a modulation of wavevector g^, in contact with a nematic solvent, and 
assume strong homeotropic boundary conditions, Sh{r±, z — 0) — — V^w(r^) , where 5n is a variation of the nematic 
director and u{rj_) is a surface displacement along its equilibrium normal z; the dimensions in the membrane plane 
are denoted by rj^. The bulk nematic solvent minimizes the Frank elastic energy Fp ~ {l/2)Kp J(fir{Wn)^ , 0|, in 
which we make the one constant approximation for the Frank constants Kp. The solution is 
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6n{v^,z)= j ^qi^^(9±)e-''^-■'■--l«-l^ (1) 

Jq±_ 

where J^^ = Jd\±/{2Tr)'^ with the hmits between an upper cutoff 27r/a and a lower cutoff 2Tr/L±, with a being 
a microscopic dimension and L± the membrane size. Substituting into the Frank energy and integrating over the 
dimension z normal to the interface, we find 

Fp = \kJ k^|3(l-e-l'^-l^^)Kq^)P « V / \u{q^)\\ (2) 

where Lz is a large distance cutoff which we take below to he — oo for an isolated membrane ot Lz — d for a 
stack of membranes spaced by d. Eq. (|^) interpolates between the ~ \<1±\'^ regime for an isolated membrane and the 
^ |(7^|*L2 long wavelength behavior in a finite system. The second expression in is an alternative approximation 
which handles properly the large and small q limits, and is much easier for calculations. We shall mostly use this form 
of Eq. (Q) in this work, since we are primarily concerned with qualitative results. To this must be added the Helfrich 
energy of the fluctuating membrane, 

Fh^^i^J d\^{Wluf + J d\^RG, (3) 

where the Gaussian curvature G integrates to zero for lamellar systems without topological defects and plays no 
further role here. 

The unusual non- analytic form of Eq. (|^) arises from the same considerations as the linear- |q| elasticity of the 
air-fiuid-solid triple line, where energy is stored in the deformation of the air- fluid surface From this energy one 
finds many respects in which membranes in nematic solvents differ from their isotropic counterparts. This behavior 
could most easily be seen by preparing mixtures of surfactant and thermotropic liquid crystal, with the latter playing 
the role of an oil (possibly with water or a cosurfactant, as is common in conventional surfactant systems, to select 
from the zoo of possible phases), and cycling through the solvent's isotropic-ncmatic transition temperature. We 
proceed by briefly describing some of these properties in order of complexity. We only consider the case of membranes 
without surface tension, where this new term is most important. In the Appendix we present the contribution due to 
the Casimir effect [Q, which is an entropic effect in the correlated fluid mediating the membrane, and show that the 
main effects of it are a renormalization of the area per surfactant head group, and the bending modulus k. 

Single-membrane properties — For these properties we take Lz = oo. The flrst obvious new effect is a qualitative 
change in the surface fluctuations. Bilayers are typically rough due to thermal fluctuations, and the combination of 
a two-dimensional surface fluctuating in three dimensions yields divergent height fluctuations. In a nematic solvent, 
however, this changes. For example, fluctuations of the surface normal are given by 



(|<5n(r)-<5n(0)r)=2fc,T / I ""'^^'^ ^ ^ log 



Kp + 2nK/a 



Kp + 27rK/r 



(4) 



>q^ KF\q± \ + nqi ttk 

where a is a microscopic cutoff. If we define the correlation length that distance along the membrane for which 
fluctuations in the normal vector n orientation are of order 1 llOfl, we find 
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C 2^K/{kBT) 

S£ ^ r (K) 

For Kp = the membrane in the isotropic non-correlated solvent is crumpled at distances larger than ^o- For Kp ^ 0, 
in a nematic solvent, increases rapidly and reaches the system size {oo) for Kpa/{2TTK) = (gW^sr _ Hence, 
for Kpg/k 3> 1, in the regime that should be identified with a 'strong nematic solvent', the membrane would not be 
crumpled at all. A typical estimate of the surfactant bilayer bending rigidity is k '--^ 5 • 10~^° J, only slightly larger 
than the thermal energy at room temperature. Taking a characteristic value for the Frank constant, Kp ^ 10^^^ J/m 
and the molecular size a ^ one obtains an estimate of order unity and, therefore, both crumpled and flat regimes 
are accessible for a membrane in a nematic solvent. 

Also of interest is the related quantity, the membrane roughness, given by the mean-square height fluctuations: 

1 



{uivY) =kBT 



+Kp\q±_\^ 

2t:k + aKp 



Kp h K log 



2t:KI 2Tr 



2ttk + L±Kp 



(6) 



where L± is the transverse membrane dimension, coming from the lower cutoff in g-space. As Kp — s- we recover, 
after expansion in powers of L±Kp, the result for a conventional membrane, (u^) ^ L'^f^, with a roughness exponent 
(■5 = 1. In a strong nematic solvent we have Cs = 1/2 and, as expected, the membrane is not as rough. 

Renormalization of bending modulus — Since the membrane in the nematic solvent is stiffer, we expect the renor- 
malization of k due to thermal fluctuations | [Tl| ] to be much reduced. There are two new sources of renormalization 
for membranes in nematic solvents: entropic, due to the Casimir effect, which we briefly discuss in the appendix; and 
energetic, due to the Berreman-de Gennes energy. Following the simple procedure outlined by Helfrich we find 



, 1 I^Qmax + Kp \ 3 L± 

1 —rr- - ^7:7 log 



(7) 



47r 

where the first correction is from the Helfrich renormalization and the second term is produced by the Casimir effect 
in the correlated solvent. In the limit KpL±/K <C 1 (isotropic solvent) the Helfrich effect returns to the usual \ogL±/a 
reduction of the bending rigidity ||ll|,|l2|. In the nematic solvent with Kp ^ it is replaced by the constant factor 
log [1 + {n/Kp a] , so that the renormalization be large or small depending on the 'strength' of nematic 

solvent. In addition, there is an increase of k, due to the Casimir effect, also logarithmically divergent with the system 
size. This result supports the intuitive expectation for the membrane to become more rigid due to the anchoring with 
the nematic solvent. 

Lamellar Phase: Helfrich Interaction — In a lamellar state we take Lz = d as the cutoff in Eq. (||), since the range 
of the solvent extends only up to neighboring membranes. We define a 'strong' nematic solvent in this context as one 
for which Kpd ^ k [which is, in fact, a much weaker condition than Kpa ^ k in Eq.(^]. To estimate this we again 
take K ~ 10"^" J , Kp ~ 10^^^ J/m, and lamellar spacings ranging from 10 — 1000 A, yielding Kpd/ n ~ 1 — 100. Since 



4 



the moduli k can be changed by adding co-surfactant and Frank constants depend on the nematic order parameter, 
it is quite easy to span the whole range from weak to strong nematic solvents. 

In stacked lamellar phases there are two well-known interactions which stabilize the lamellar phase: electrostatic 
stabilization and stcric interaction. The electrostatic effect yields an interaction energy per unit area of F/A ^ 
kgT / [LBd)^ where is the Bjerrum length and d is the membrane separation. We shall not consider this interaction, 
exploring instead the more interesting statistical effects of electrostatically screened membranes In an isotropic 
solvent these effects lead to the steric stabilization [^, giving F/A ^ [k bT)'^ / [nd^) . This steric interaction arises 
due to the divergence of the height fluctuations of a single membrane, ^ L?f^; membranes in a stack experience 

collisions with a characteristic length between collisions governed by Qs- Since is positive we expect a steric 
interaction in the nematic solvent as well, but with a different character than for the standard Helfrich interaction. 

Rather than going through a calculation similar to Helfrich's, we content ourselves here with a scaling-type analysis 
to obtain the d-dependence of the steric stabilization of membranes in a nematic solvent. We first calculate the height 
fluctuations, using the equipartition theorem for the membrane Hamiltonian, given by Eqs. (|) with L^ = d and (|), 
and integrating over wave vectors qj_ 



/ i fcj^r l + [K^d/K]{l + 2qod) 

^"^"■^ ^ ^ 47rA^g2 (l + [K,d/K])^ ' 



where go = 2tt/Lj_ is the low-g cutoff in the membrane plane (we have ignored subdominant terms logarithmic in qq). 
Now, the membrane will be sterically stabilized when the height fluctuations are of order the layer spacing, (u^) ~ d^. 



This determines qq. For a weak nematic solvent, Kpd <C k, we recover the Helfrich result qod « y^hgT/ k. For a 
strong nematic solvent we find 

k T 

M)' = T^(l + 29od); (9) 
AnKf-d 

Since, typically, k > kgT we have Kpd ^ kgT, in which case w {kuT / A-kK pY^'^ d~^^'^ . This defines a new length 
^ — q^^ ^ the in-plane correlation length for height fluctuations, or mean distance between collisions. Now we may 
compute the interaction energy. Crudely, the pressure P due to undulations may be calculated as that of a gas of 
sterically interacting 'discs' of dimension ^. This yields P ^ kp,T / [d^"^). Using our result for ^, and realizing that the 
free energy per unit area is then given by F/A ~ Pd, we arrive at 

^ = (10) 

which should be compared with the d^^ and d^^ behavior of, respectively, the standard Helfrich and electrostatic 
stabilizations. 

The conclusion is that, for electrostatically screened membranes in a nematic solvent, much closer lamellar packing 
can be achieved, which again confirms the intuitive expectation for more rigid and flat membranes. The resulting 
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"smectic" lamellar phase will have a different layer compression modulus of the corresponding Landau-Peierls elastic 
energy 



(11) 



where K = n/d is the layer bending modulus in the absence of nematic solvent. In writing this we have ignored the 
\q±\d term in the denominator of Eq. (|^), since smectic elasticity is concerned with wavelengths much larger than the 
smectic spacing, \q±\d ^ 1. The compression modulus is B = B — C^X: where B is the bare compression modulus, 
and B includes the renormalization due to the coupling between solvent composition and layer spacing Jl6t . B is 
given essentially by the pressure of the gas of colliding membranes, which scales as B 1/rf'', where p = 2, 3,4 for 
electrostatic, standard Helfrich, and nematic-solvent membranes Eq. ([lo|) . The compression modulus may be measured 
by, for example, small angle scattering [ p^ . 

More useful information can be extracted from the line shape of the diffusion scattering peak, following from 



the Landau-Peierls energy (11). As in ordinary smectics, the structure factor behaves as, for example, S(0,qz) ~ 
{Qz — Id)^"^^^^ with the usual Caille exponent |p7|| , 

^-gg -TH- (12) 

Because Kp ^ if in strong nematic solvents, 77 should strongly decrease when the solvent undergoes a nematic 
transition, leading to a more rapid decay of the structure factor, as noted above. Notice also that the exponent 
depends on B, which changes its qualitative dependence on d as one moves into the strong solvent regime. By 



measuring 77 and the penetration depth A = ^JB /{K + Kp) (which may be extracted by the corrections to the low-g 
behavior of iS'(0, qz) p^) one may determine both B and K + Kf by systematic dilution and temperature experiments. 

Unbinding Transition — Upon lowering the temperature into the nematic phase, the steric repulsion energy F/A at 
fixed d drops by a factor /«/// ~ K/{dKF) <C 1. This dramatic decrease should affect the unbinding of layers |l9| , p0| . 
Within a Flory-type theory, Milner and Roux showed that one can write the free energy per volume / of a stack of 
bilayers as |l^] 

(") 

where c is a numerical constant, 6 is the bilayer thickness, and x accounts for contributions to the second virial 
coefficient from other than steric {i.e. typically van der Waals) interactions. Here <j) 5/dis the surfactant volume 
fraction, and the (jr' term follows from the isotropic-solvent steric interaction. When the solvent undergoes a transition 
into the nematic phase the (j)^ term should be replaced by 
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where c' is another numerical constant. Since Kpd can be much larger than k, the result is a smaller repulsion and a 
smaller preferred interlayer spacing. 

At fixed (f) the characteristics of the unbinding transition are: at low x ^ single bound phase exists, while Xc = 
icksTcj)/ (k5'^) marks the spinodal line at which the system phase separates into bound (0 7^ 0) and unbound (</> = 0) 
phases. 

For nematic solvents the spinodal line is given by 

which allows for the possibility of, for example, phase separation by quenching the solvent into a nematic state. The 
resulting dynamics would be very complicated, due to the simultaneous nematic coarsening and phase separation. 

Conclusions — In summary, a membrane in a nematic solvent should be much stiffer than in an isotropic solvent, 
leading to its different scaling behavior. Layered systems are in this stiff regime when Kpd/n ^ 1, which should be 
experimentally realizable. The stiffening can be seen in several quantities, such as the correlation length ^0 for surface 
normal fluctuations; the roughness exponent ^5; and intermembrane interactions, in which the standard Helfrich 
interaction changes its dependence on the intermembrane spacing d. In addition to the effect of the bulk nematic 
elastic energy, a membrane is also affected by the entropic Casimir-like effect of fluctuations of the coupled director 
field. We have not yet considered the case of smectic solvents, but it is straightforward to show, by arguments very 
similar to those used in deriving Eq. (Q), that there is a simple renormalization of the membrane elasticity modulus, 
Ka = K + 2\/KsBsd, where the subscript s refers to the solvent smectic elastic constants; and the Helfrich interaction 
is the same as in the isotropic case, with k replaced by Kr. 

Interesting effects are expected in the presence of an external magnetic field H, which provides a "mass" for the 
nematic director fluctuations. While the effect of magnetic field on a membrane is, in principle, the same as for an 
isotropic solvent, the anisotropy of the diamagnetic susceptibility of a bilayer membrane should be negligibly small 
compared to that of a bulk nematic liquid crystal. It is straightforward to show that the application of a magnetic 
field H along the layer normal yields a term in the free energy Eq.(|2|) (for a single membrane in an infinite system) 
proportional to q\\J q\ + ^h^, where = {^Kp / (xaH'^)Y^'^ is the standard magnetic coherence length. This term 
further reduces the height fluctuations of the membrane to ~ logi^/a, which leads to a very weak steric 

Helfrich repulsion F/A ^ l/d^, and suppresses the Landau-Pcicrls instability in favor of Bragg peaks at the smectic 
wavevector. 

It seems fairly straightforward to perform experimental checks on the described system, by mixing a thermotropic 
nematic with a small concentration of surfactant, choosing its hydrophobic part to be closely related to mesogenic 
molecules. Addition of a small amount of water would further stabilize the bilayer membrane structure. All our 
arguments suggest that it would be very difflcult to create curved micellar structures in the nematic solvent (spherical 
micelles, for example, would have to create a topological defect in the nematic field around them, due to the radial 
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director anchoring). Instead, we expect the formation of rather flat bilayers even at very low concentrations and dense 
lamellar and sponge phases with more surfactant/ water, with the morphology driven by the elastic energy effects in 
the mediating nematic solvent. 

We appreciate useful discussions with M. Warner and the support and practical advice from the Polymers & Colloids 
group of Cavendish Laboratory, which allowed us to observe some of the described effects in practice. This research 
has been financially supported by Unilever-PLC (EMT) and the EPSRC (PDO). 
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Appendix: Casimir Effect — In addition to the energy stored in the director field there is an entropic contribution to 
the membrane free energy due to the analog of the Casimir effect, calculated for liquid crystals by Ajdari, et al. j|] 
and Li and Kardar From the results of Li and Kardar, the entropic contribution to the free energy per unit area 
of a membrane fluctuating above a flat surface a distance d away is [Q 
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1 + 4.C,(1)Q^ 



log 



47rCi(2) 



where oi = 0.04792, and Ci(C) is given by Eq. (2.19) of Ref. [4b] and, generally, are very small. To obtain this we 
have expanded Li and Kardar's results, which hold for an arbitrary surface, in a gradient expansion in the membrane 
fluctuation u{rj_). The term multiplied by ai is essentially the classic Casimir attraction between the plates, and 
contributes a small renormalization to the existing attractive interactions, which are typically van der Waals [ p2[ . 
The term renormalizes the surface tension and, since a surfactant system in solution adjusts its area per 

head group S to retain equilibrium and satisfy vanishing surface tension po| , leads only to a slight decrease in S. The 
last term, written in the single- membrane limit d — > oo, renormalizes the membrane bending rigidity k by a relevant 
logarithmic term (see Eq. (|^), where we neglected the small correction Ci(2) ^ 10"'^), which should be compared with 
the value [kBT/AT:] log{L±/a) found for the renormalization due to thermal fluctuations of membranes in isotropic 
solvents lllll. 
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FIG. 1 Director field n near a fluctuating membrane, which imposes homeotropic anchoring 
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